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ABSTRACT 

We use high-resolution N-body simulations, combined with a halo occupation model of galax y bia s , to in - 
vestigate voids in the galaxy distribution. Our goal is to address the "void phenomenon" of PeeblesI (1200 ll) . 
which presents the observed dearth of faint galaxies in voids as a challenge to the current cosmology. In our 
model, galaxy luminosity is determined only as a function of dark matter halo mass. With this simple assump- 
tion, we demonstrate that large, empty voids of ~ 15 /i~' Mpc in diameter are expected even for galaxies seven 
magnitudes fainter than L,. The predictions of our model are in excellent agreement with several statistical 
measures; (/) the luminosity function of galaxies in underdense regions, (ii) nearest neighbor statistics of dwarf 
galaxies, (///) the void probability function of faint galaxies. In the transition between filaments and voids in 
the dark matter, the halo mass function changes abruptly, causing the maximum galaxy luminosity to decrease 
by ^ 5 magnitudes over a range of ~ 1 /;"' Mpc. Thus the boundary between filaments and voids in the galaxy 
distribution is nearly as sharp for dwarfs as for ^ L, objects. These results support a picture in which galaxy 
formation is driven predominantly by the mass of the host dark matter halo, and is nearly independent of the 
larger-scale halo environment. Further, they demonstrate that ACDM, combined with a straightforward bias 
model, naturally predicts the existence of the void phenomenon. 
Subject headings: galaxies; halos — large scale structure of the universe 



1. INTRODUCTION 

Galaxy redshift surveys have revealed a complex network 
of clusters, filaments and walls. They have also demon- 
strated that the expansive regi ons between these structures 
are nearly bereft of ga l axies (seelGregorv & ThompsonlI978l; 
iKirshner et ap [198It_ l^elev et al.1 Il99j for early work, 
and 'Hovle & Vo2elev"2004'; 'Croton et al."2004t IConroy et al.) 
|2005; Patiri et al. 2006a; Tinker et al. 2007a for results from 
recent large-scale surveys). The surprisingly large size of 
voids — up to ~ 30 h~^ Mpc in diameter — and their apparent 
depth, in terms of their luminosity density, have begged ques- 
tions about their formation mechanism. Why do faint galax- 
ies avoid the voids defined by their brighter brethren? Further, 
the few objects found in and near underdense regions appear 
to repres ent a fair sample of the overall galaxy population. 
Peebles (1200 U hereafter POl) presented these observations 
as "the void phenomenon". The straightforward argument 
is that the deepest voids in the dark matter distribution pre- 
dicted by standard ACDM are never completely empty; they 
still contain mass and significant numbers of low-mass halos. 
Therefore, voids should be the preferable environment of faint 
dwarf galaxies. However, POl demonstrates that dwarf galax- 
ies tend to be distributed similarly to bri ghter spirals, and the 
local void exhibits a paucity of dwarfs. Peebles (2007) further 
shows that the edge of the local void has the same sharp edge 
for both bright and faint galaxies. If halos exist in voids, then 
where are the galaxies? 

In this paper we present a quantitative model of galaxy 
bias that allows us to directly address this question. 
We use the Halo Occ upation Distribution (HOD; e.g . 
Seliak 2000; Scoccimarro et al.ll200lHCooray & Shethll2002t 
Berlind & Weinbera r2002h to specify the relationship be- 
tween galaxies and dark matter. In the HOD framework, 
this relationship is quantified by P{N\M), the probability 
that a halo of mass M contains A^ galaxies of a given sam- 
ple. Galaxy samples can be defined by luminosity, color. 



star formation rate, or any other galaxy property. For each 
sample, P(N\M) will be different but can be calibrated with 
observational data. The HOD has emerged as the dom- 
inant tool for interpreting clustering measurements for a 
wide range of redshifts and galaxy cla s ses (see, for a tip 
of the iceberg, jZ ehavi et al. 2004, 2005; Zheng et al. 2007; 



van den Bosch et al...2 007 ; Tin ker et al...2007b.a; .Chen.2007 ; 
White et al.ll2007l: iPad manabh an et al.ll2008l) . The results ob- 



tained from analysis of observational data are in excellent 
agreement with the results of semi-analytic models, hydro- 
dynamic cosmological simulations, and high- resolution coUi- 
sionless simiilation s (.Kravtsov et al...2004: Zheng et alJl2005t 
IConrov et al.ll2006l) . Our knowledge of the mapping between 
galaxies and halos is now well-established and can be confi- 
dently extended into regimes in which observational data are 
lacking. Here we take HOD results calibrated on observed 
galaxy samples that occupy > 10" /i~' M0 halos and extrap- 
olate to lower masses in order to model the distribution of 
dwarf galaxies. We will demonstrate that the predictions of 
this model are in excellent agreement with the observations of 
POl: deep, wide voids are expected within hierarchical struc- 
ture formation, even for the faintest of galaxies. 

One of the difficulties in directly addressing the void phe- 
nomenon is that the problem has not been clearly defined. 
The definition of a dwarf galaxy varies in the literature, and 
the predictions of ACDM in this regime have not been ro- 
bustly specified. To elucidate the problem, POl compared 
the spatial distribution of "ordinary" galaxies to samples 
of "test" d warf galaxies within the Optical Redshift Sur- 
vey (ORS; [Santiago et al.l Il995h . The number density of 
the ordinary sample is comparable to the number density of 
Mr < -16 galaxies ' in the Sloan Digital Sky Survey (SDSS; 
iBlanton et al.ll2005h . Thus to give our model a fixed target we 

' For brevity, all galaxy magnitudes assume a Hubble constant h = 
Ho /\00= 1. All r-band magnitudes have been ir-corrected to z = 0, as with 
the lBlanton et all 120051) luminosity function. 



stipulate that the void phenomenon begins at galaxies fainter 
than Mr = -16, roughly Small Magellanic Cloud-type galaxies 

and fainter. 

In [Tinker et alj ( l2007al) we used the HOD to make pre- 
dictions for void probability statistics for galaxies brighter 
than M, = -19. Our predictions were in excellent agree- 
ment with our measurements f rom SDSS Data Release 4 
jAdelman-McCarthv et alj l2006l) . Thus we concluded that 
there is no conflict between the observed voids and those 
predicted by ACDM for galaxies down to this magnitude 
(roug hly 0.2L,, given an M» of -20.44 from Bl anton et alj 



jgh 

m. 



l2003h . M ore importan t ly for t he model presented herein, the 
results of ' Tinker et al.l (l2007ah support a model in which the 
luminosities of galaxies are determined entirely by the mass 
of the host halo, independent of the environment in which the 
halo formed. This allows us to robustly predict the distribu- 
tion of void galaxies within a model that connects galaxies 
with halos through the halo mass only. 

Several other studies have investigated galaxy voids 
and void gala x ies th rough a variety o f tec hniques. 
iMathis & Whit3 (|2002|) and iBenson et all (l2003h . using 
semi-analytic models, concluded that galaxies of M^ ~ -18.5 
avoid the voids of the brighter galaxies. These studies 
were limited by numerical resolution and were not able to 
model the dwarf galaxies on which POl based his argument. 
[Furlanetto & Piran (2006) constructed an analytic framework 
within which to explore the dependence of voids on galaxy 
luminosity, bui lding on earlier work modeling voids in the 
dark matter by ISheth & van de Weygaerj (l2004l) . However, 
no tests with numerical simulations or comparison to obser- 
vational results were made. The model of Furlanetto & Piran 
(J2006) is an excellent tool for understanding the trends seen 
in the data, but detailed analysis of voids requires a simulation 
to produce the correct halo distribution. As we will show, the 
change in the halo mass function along the edge of a halo 
is of critical importance in the resulting void distribution. 
Additionally, observational data are all measured in redshift 
space, which is not incorporated into these analytic models. 

This paper is organized as follows: In §2 we present our 
HOD model and the numerical simulations employed to cre- 
ate mock galaxy distributions. In §3 we present the predicted 
void distribution, making quantitative comparisons to several 
sets of observational data: (a) the luminosity function of void 
galaxies, which h as been measured down to Mr = -14 by 
iHovle et all (l2005h . (b) the nearest neighbor statistics of POl, 
which probe galaxies fainter than M,. = -16, and (c) the void 
probability function of galaxies down to Mr = -17, for which 
we present new measurements from SDSS Data Release 6 
dAdelman-McCarthv et al.|[2007h . In §4 we discuss our re- 
sults. In all calculations, we assume a flat, ACDM universe 
with (a«,cr8,/j,n, fife) = (0.3, 0.9, 0.7, 1.0,0.04). 



2. METHODS 



2.1. 



Halo Occupation Models 

Halo occupation models generally contain two important 
mass scales. The first is the mass at which halos become suf- 
ficiently large to host a single galaxy brighter than a specified 
luminosity threshold. The second is the mass at which ha- 
los become massive enough to host additional satellite galax- 
ies. Results from both theory and observation have demon- 
strated that these masses scale in a straightforward fashion 
with galaxy luminosity, allowing for confident extrapolation 
to the dwarf galaxies on which we focus. 



TABLE 1 
Minimum Halo Mass as a Function of 
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Note. — Galaxy densities are in units of 
r' Mpc)-^ . All masses are in units of ft-' Mq . 



Functionally, the HOD breaks down the occupation of 
galaxies within halos into two distinct components: central 
galaxies and satellite galaxies. We use the standard parame- 
terization for central galaxies brighter than a given magnitude 
limit: 
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where M^in is the mass at which a halo has a 50% probabil- 
ity of having a central galaxy above the defined luminosity 
threshold, and CTiogM is the width of the transition between 
and 1 galaxies, physically representing the scatter of mass at 
fixed luminosity. The number of satellite galaxies scales as a 
power of the host mass, with a cutoff scale set by the central 
occupation function. 
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(A^sat)M = (A^cen)M X ( — 
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where Mi is the mass scale at which a halo has, on aver- 
age, one satellite brighter than the defined magnitude limit. 
The inclusion of (A^cen)M in equation dU ensures that the cen- 
tral galaxy is the brightest galaxy within a halo. The forms 
of these equations are in good agreement with the resu lts of 
semi-analytic models of galaxy formation dZheng et alj20 05). 
high-resoluti on coUisionle ss N-body simulati ons that re solve 
substructure (iKravtsov et a l. 2004; Conrov e taP l2006h . hy- 
drodynamic simulations (Zheng et al. 2005), and analytic 
models of halo substructure (Zentner et al. 2005). 

It has been demonstrated that the parameters of the occupa- 
tion functions are nearly self-similar with luminosity, i.e., that 
ttsat ~ 1 and Ml /M^nin ~ 20. These results have been found 
both in the theoretical results listed above and analyses of 
clustering measurements from numerous surveys at multiple 
redshifts (IZehavi et a l. 2005; Zheng et al. 20()7t iTinker et all 
l2007bt Ivan den BosTh et al...2007.) . High-resolution N-body 
simulations have demonstrated that the subhalo mass func- 
tion is nearly self-similar with parent halo mass, in good 
agreement with the HOD results obtained from observation s 
(IKravtsov et al.ll2004l; iGao et alj|2004t iDe Lucia et alll2004l) . 
These results extend down to subhalos of M ~ 10^ /j~' M© 
(iDiemand et al.ll2007l) . With the now well-established con- 
nection between s atellite galaxies and substructure within 
dark matter halos dKravtsov et al.ll2004t IConroy et aLJ 120061: 
IWeinberg et al.l l2006). the self-similarity of the HOD with lu- 
minosity is expected to extend far down the mass function to 
halos that host dwarf galaxies. In addition, because we are 



TABLE 2 
Properties of the Simulation Set 
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Note. — Lbox is in units of h^' Mpc. e 
is the Plummet foice softening length. Np is 
the total number of particles in the simula- 
tion. Magnitudes in column 5 represent the 
lowest magnitude bin modeled by each simu- 
lation. 



interested in voids, the details of the (A^sat)M are largely iiTel- 
evant; as long as the fraction of galaxies that are satellites is 
roughly correct, t he distribution of voids is not affected b y the 
details of (Mat)M (IConrov et al.ll2005l: [Tinker et al.ll2006h . 

For each luminosity threshold, we set the val ue of Mmn by 
matching the space density of galaxies from the lBlanton etal.1 
( l2005h luminosity function, ie. 



/•OO 7 

'i>{Mr)dMr = / {{N,,„)m+ (Mat)M) ^ dM , (3) 
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where dn/dM is th e hal o mass function, for which we use 
[Tinker et alj ( |2008|) . The lBlanton et all (l2005h result is accu- 
rate down to Mr < -12 galaxies. We set M\ = 22Mj„i„ and 
ttsat = 1 for all samples. For criogM we choose a value of 0.15, 
which is consistent with the results of Tinker et al. (20073) 
for Mr < - 19 samples. This is somew hat smaller than that 
obtained bv lvan den Bosch et alj (l2007h . but we note that the 
scatter between mass and luminosity is expected to be larger 
for blue-selected galaxy samples (which Ivan den Bosch et al.l 
I2OO7 analyze). Tinker et al. (2006) concluded that the values 
of CTiogM ranging from to 0.5 produced little effect on the 
void distribution. We will test the dependence of our results 
on (TiogM in §3. Table 1 lists the value of Mmin for each lumi- 
nosity threshold. 

2.2. Numerical Simulations 

To investigate the distribution of voids and void galaxies 
with our HOD models, we populate a series of high-resolution 
N-body simulations of various box sizes kindly provided my 
M. Warren. The parameters of the simulations are listed in 
Table 2, and we will refer to each simulation by its box size 
in /z"'Mpc: L96, L192, and L384. The simulations allow 
us to probe halos robustly down to ~ 10'" /i~' Mq, the mini- 
mum mass scale of -12 galaxies. These sim ulations were per- 
forme d using the hashed oct-t ree code of . Warren & SalmonI 
(Il993h . and were described in IWarren et al. j 20061). We use 
the spherical overdensity halo finder of iTinker et al.l (l2008l) 
to identify halos in the simulation, defining halos with an 
overdensity A = 200 with respect to the background matter 
density. The mass resolution of the highest resolution simu- 
lation is 18 times hig her than the Millennium Simulation of 
ISpringel et al.l (l2005b and 75 times higher than the simula- 
tion coupled to the sem i-analytic galaxy formation model of 
lMathis&Whita(l2002h . 

To populate each simulation with galaxies, we use a Monte 
Carlo approach based on equations ([U and (|2|i, assuming a 



nearest integer distribution for central galaxies and a Pois- 
son distribution of satellite galaxies. The nearest integer ap- 
proach is standard for central galaxies due to the limit of one 
central galaxy per halo. Poisson statistics provide an excel- 
lent description of the distribution of satellites in simulations 
( Kravtsov et al.ll2004t IZheng et anl2005) and in observations 
dLinetal .1120041) . By considering the HOD parameters as a 
function of luminosity, the full conditional luminosity func- 
tion of each halo can be easily determined. The number of 

centrals in a binM, to M, + AM, is (A^cen)M ''- (A^cen)M '"''^*''*- 
The number of satellite galaxies in a magnitude bin is calcu- 
lated in the same manner 

Central galaxies are placed at the center of the halo, and 
satellite galaxies are distributed randomly with the density 
profile of N avarro et al. (1997), with the concentration-mass 
relation of Bullock et al. | (,2001.) (using the u pdated parameters 
of the model listed in IWechsler et al.ll2006l) . Satellite veloci- 
ties, relative to the halo motion, are taken from a Gaussian dis- 
tribution with 1 -dimensional dispersion of cr^ = GM2(K)/2/?200- 
We populate each simulation down to a limit of 100 particles 
per halo. The 384 /z"' Mpc simulation is ideal for M,- < -20 
samples, the 192 /i"' Mpc sample resolves M, < -17 samples, 
while the 96 /i~' Mpc simulation probes galaxies as faint as 
Mr = -14. Although Mmin for Mr = -12 is 127 particles, we 
must resolve halos below M„^i„ to fully account for the scatter 
between mass and luminosity. 

3. RESULTS 

3.1. The Cosmic Web in the Faintest Galaxies 

Figure [Upresents a slice through L96 that is typical of the 
structure in these simulations. The galaxy positions are in 
redshift space, assuming the distant-observer approximation 
with the y-axis being the line of sight. Along the filaments, 
L* and dwarf galaxies cohabitate, along with all luminosi- 
ties in between. The transition between filaments and voids 
is very sharp, however. Deep voids of ~ 10 /!~'Mpc in di- 
ameter, empty of even the faintest resolved galaxies, are not 
uncommon in this simulation. The deepest void, centered at 
{x,y) = (80,48) /i"' Mpc, is nearly 18 /1"' Mpc along its longest 
axis. The transition from the main filament in the slice (run- 
ning along the y-axis) into the void on either side is sharp for 
both bright and faint galaxies, resem bling the galaxy distribu- 
tion around the local void shown in iPeeblesI (|2007|) . Around 



the local void, between 5 to 7 /;"' Mpc from the void center, 
the number of galaxies Mb < -18 and Mb > -18 drops es- 
sentially to zero. The abrupt transition from filament to void 
is ubiquitous in the cosmic web, even for galaxies hosted by 
10'" /z~' M0 halos. We note that halo velocities are usually 
directed along the axis of the filament, thus choosing the x- 
axis as the line of sight does not noticeable change the galaxy 
distribution plotted in Figure [T| 

3.2. The Luminosity Function of Void Galaxies 

The most straightforward quantifiable test of our model is 
the number of galaxies located withi n voids, and their d istri- 
bution as a function of luminosity. iHovle et al] (l2005l) cre- 
ated a sample of "void" galaxies in the SDSS by identify- 
ing galaxies with local densities below a critical value of 
5 = 5p/p = -0.6 in a top-hat sphere of 7 /;"' Mpc. Densities 
were calculated with respect to galaxies with M, < -20.5. We 
have constructed a sample of void galaxies in our simulations 
following these same criteria. All densities are calculated in 
redshift space assuming the distant observer approximation. 
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Fig. 1 . — Slice through the galaxy population created from the halos within the 96 ft"' Mpc simulation. The depth of the slice is 6 ft*' Mpc. The point size 
scales with the luminosity of the galaxies, r-band magnitudes are as follows: blue = — 14and— 15; green = —16 and— 17, red = —18 and— 19, orange are all galaxies 
—20 and brighter. 



Figure [2] com pares the vo i d lum inosity function of our HOD 
model to the iHovle et aH (12005.) measurements. Our model 



is in excellent agreement with the data, not only reproducing 
the overall abundance of void galaxies but also the decrease 
in the value of M, in the void luminosity function relative to 
the overall luminosity function. 

Interestingly, we find some void galaxies with magnitudes 
as bright as M,- = -21 .5, just as measured in Hoyle et al. These 
objects are not scattered into void regions due to redshift 
space distortions, but are in intrinsically low-density regions 
as defined by the M,- < -20.5 galaxies. The minimum mass 
scale for Mr < -21.5 galaxies is ^^ 1.8 x lO'-' /!"'M0, which 
are never found in 5 < -0.6 regions in the dark matter dis- 
tribution. When using the dark matter particles to obtain the 
local density around each galaxy and recalculate the void lu- 
minosity function, there are no void galaxies this bright. Thus, 
bright void galaxies are due to stochastic biasing of L* galax- 
ies, both from Poisson fluctuations in the number of halos and 
fluctuations in the number of L > L, galaxies per halo. This 
creates a few regions in which the dark matter density is above 
the density threshold but the galaxy density is below it. 

3.3. Nearest Neighbor Statistics 

POl used nearest neighbor statistics to probe the relative 
distribution of "ordinary" galaxies (Mr < -16) to dwarf galax- 
ies (Mr > -16). If dwarf galaxies preferentially occupy void 



regions relative to their brighter counterparts, their distribu- 
tion of nearest neighbors will show a significant tail out to 
large neighbor distances Rn^. To circumvent the problem of 
different galaxy samples having different mean space densi- 
ties, for each test object the nearest neighbor in the control 
sample is found. This distribution is compared to the distri- 
bution of R„b of the control sample to itself. POl found that 
the distribution of nearest neighbors for test and control ob- 
jects are essentially the same, indicating that dwarfs avoid the 
voids defined by the ordinary objects. 

In Figure [3] we present several examples of the cumula- 
tive distributions of Rnb taken from L96. Here we have set 
our control sample to be Mr = [-16,-18), and our test sam- 
ple are galaxies with Mr = [-14,-16). The control and test 
samples are substantially different in terms of their space den- 
sities, but are not too dissimilar in the halo masses that they 
probe. To test any systematics due to these choices we have 
an additional sample of "bright" objects containing all galax- 
ies Mr < -18. The largest sample in POl is roughly 1000 
(/i"' Mpc)"*, so we have broken our mock galaxy distribution 
into 9"^ equal-volume cubes of 1225 (/i~' Mpc)-' to test for cos- 
mic variance in this statistic. The bottom right panel shows 
the results from the full simulation. We define R,c as the dis- 
tance to the nearest control galaxy for each test galaxy, while 
Rbc is the distance to the nearest control galaxy for each bright 
galaxy, and Ra- is the distance for each control galaxy to the 
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Fig. 2. — Comparison betw een the luminosity function of void galaxies 
measured bv lHovle et alj 120051) and predicted by our HOD model. The HOD 
measurement, shown with the open squares, is a weighted average of results 
from the 96 /i"' Mpc and 192 ft"' Mpc simulations (noting that all points at 
Mr > —17 are from the 96 h^^ Mpc simulation). En'ors are estimated by jack- 
knife sampling of each simulation into octants. The solid line is the Hoyle 
et al. data, with the shaded region representing their error bars. For compar- 
ison, the full luminosity function of all galaxies from lBlanton et all f2003) is 
shown with the dotted curve. 



nearest galaxy within the control sample. 

The shaded region in each panel of Figure |3] approximates 
the results from the largest dwarf sample from POl (his Fig- 
ure 4, top set of curves). The upper limit of the shaded region 
is set by /(< Rtc) while the lower edge is /(< Rcc)- POl de- 
termined the ratio of mean /?nb values to be {Ric)/{Rcc) = l-l- 
This is nearly identical to the overall value of 1.06 obtained 
from the L96 simulation. The scatter in this statistic is quite 
large, with a variance in mean distance ratios over the 9^ sub- 
samples of 0.23. In many examples in Figure |3] /(< Rcc) 
extends to distances significantly beyond /(< R,c). Compar- 
ing the R,c distributions to the nearest neighbor statistics for 
the bright objects to the test sample, Rtc yields similar re- 
sults. The mean distance ratio is 1.19 with a scatter of 0.55, 
owing to the lower number density of brighter objects. While 
the agreement between model and data are encouraging, it is 
likely fortuitous given the large cosmic scatter in the POl mea- 
surements. However, we conclude that the prediction of the 
model supports the picture that dwarf galaxies and brighter 
galaxies have similar spatial distributions in low-density re- 
gions. The slight difference in nearest neighbor distributions 
between dwarf and regular galaxies i s also i n agreement with 
the observational results of iLee et al] (l2000h . 

3.4. Void Probabilities as a Function of Luminosity 

The void probability function (VPF, denoted PQ{r)) is 
defined as the probability that a rand omly placed sphere 
of radius r contains no galaxies. In Tin ker et aH (12007 ah 
we measured the VPF for galaxies as faint as Mr = -19 
from Data Release Four of the Sloan Digital Sky Survey 
dAdelman-McCarthv et al.l [20061) . Fainter samples were not 



used because cosmic variance errors become large and galaxy 
clustering measurements are not available at Mr > -18. In 
Figure|4]we present new measurements of the VPF from DR6. 
The measurements are in magnitude bins, 1 magnitude wide 
(referenced by their lowe r limit). We follow the procedures 
outlined in detail in Tinker et al.l (12007 a) for both measure- 
ments and for comparing the HOD predictions to the data. 
We use the HOD mocks to es t imate e rror bars on the data, 
also discussed in iTinker et al.l (l2007ah . The increased sky 
coverage of DR6 attenuates (but does not eliminate) the cos- 
mic variance considerations for fainter samples. Though we 
are not able to probe void statistics for dwarf galaxies (in our 
definition), we are able to make VPF measurements down to 
Mr = -17, a two magnitude improvement on our previous re- 
sults. We use these data to test the robustn ess of our HOD 
model and to extend the conclusions of Tinker et al.l (l2007ah 
discussed in §1 to lower luminosities. 

Rather than presenting VPFs for each luminosity bin, we 
consolidate the results by presenting void probabilities at a 
fixed distance scale as a function of luminosity. The open 
circles in Figure |4] show Pair) at the mean intergalactic sep- 

--1/3 

aration for each magnitude bin, r = Ig = ng . Squares and 
triangles represent Po{r) at r = 1.3/i, and r = l.6lg, respec- 
tively. At luminosities below L^,, the void probability at 
any multiple of Ig is essentially independent of luminosity. 
Brighter galaxies have a lower probability of finding a void, 
thus, in the scaled distance r/lg, the brightest galaxies have 
the smallest voids. This has al so been seen in the 2dFGRS 
data ([yon Benda-Bec kmann & M iiller 2008). This is partially 
a result of the fact that brighter galaxies have fewer satel- 
lite galaxies (a higher satellite fraction results in an increase 
in Mmin in order to match the number density of the galaxy 
sample), but it is mainly due to the fact that Ig for samples 
on the exponential tail of the of the luminosity function in- 
creases much more rapidly than the bias of those galaxies. 
For Mr < -22 galaxies, Ig = 33 /i"' Mpc. 

The HOD model predictions are shown in Figure |4] for 
all three simulations. The difference in the amplitude of 
P()(r) between L96 and the larger simulations is consistent 
with cosmic variance due to the small volume. Although the 
model is not calibrated on clustering measurements as done in 
[Tinker et alJ (l2007ah . it matches the observed void probabili- 
ties as a function of luminosity and scale. The constraints on 
Mmin are driven primarily by fig rather than two-point cluster- 
ing, therefor e Mmin and the resulting VPF are similar between 
[Tinker et al; (2007a) and the model presented here. 

Although we cannot probe the VPF for galaxies as faint as 
M,- = -14 the success of our model in matching Poir) down to 
M, = -17, the void luminosity function down to Mr = -14.5, 
and nearest neighbor statistics down to M,. = -14 strongly sug- 
gest that our extrapolation of the HOD below halo masses of 
10" h~^MQ is robust, and the resulting galaxy distribution 
represents a complete picture of the structure within under- 
dense regions. 

3.5. The galaxy structure within voids 

In the VPF results of Figure H] Po(l-6/g) ^ constant for 
-17 > Mr > -20 galaxies in both the model and the data. In 
this magnitude range, voids are self-similar At fainter mag- 
nitudes, the void probability monotonically increases with 
decreasing brightness. This implies that voids are not self- 
similar for dwarf galaxies, and that the structure of the comic 
web itself plays a role in the distribution of dwarfs in voids 
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Fig. 3. — Nearest neighbor statistics for "test" galaxies, Mr = [—14,-16) (solid histograms), "control" galaxies, M, = [-16,-18), (dashed histograms), and 
"bright" galaxies, Mr < —18, (dotted histograms). For each galaxy in a sample, R^b is the distance to the nearest control galaxy in redshift space. The first eight 
panels represent random subvolumes of the 96 /r' Mpc simulation that match the volume of the sample in POl. In each panel, the number in the bottom left is 
the ratio of mean R„(, values, {R,c)/{Rcc)- The number in the top right is {Ric)/{Rbc)- The bottom right panel are the overall statistics for the entire box. The 
shaded region represents the observational results from POl. 



(as seen in Figure[TJ. 

The results of Figure [T] and the trend in Poir) for Mr > -17 
galaxies can be interpreted through the change in the halo 
mass function within the void itself. Figure |5] plots the max- 
imum halo mass as a function of distance from the center of 
the three deepest voids in L96 (all with radii ~ 10 h~^Mpc, 
and central densities 6 < -0.9). At 10 /i~' Mpc from the cen- 
ter of the void, halos between lO'^-lO'^ /z~' Mq can be found. 
These halos house the L, galaxies that 'define' the edge of the 
void. But between 10 /;"' Mpc and 5 /;"' Mpc from the void 
center, the maximum halo mass drops nearly three orders of 
magnitude. At the very centers of the voids, only M ~ lO' 
/i~' M© halos are present. The points in Figure |5] show Mmn 
as a function of luminosity from Mr < -21 to M,- < -10 sam- 
ples^. The points are placed along the curve to demonstrate 
how large the void would be in each luminosity. The steep 
drop in the maximum mass at 7? ~ 6 /i~' Mpc causes the voids 



^ We caution that the lBlanton et all 1200511 luminosity function is calibrated 
only down to —12, so results at lower luminosities are an extrapolation. 



to have sharp boundaries in the galaxy distribution over a wide 
range of galaxy luminosity. Figure |5] implies that between 5 
and 6 /?"' Mpc from the void center, the brightest galaxy lu- 
minosity will plummet by nearly 5 magnitudes. The galaxies 
in the inner half of the void (if there are any) would be well 
below the magnitude limit of the ORS data with which POl 
calculate their statistics. .Peebles! (|2007|) finds that the abrupt 
transition into the local void, inside of which no galaxies are 
observed, occurs at /? '--^ 5-7 /i"' Mpc from the void center He 
further demonstrates that this boundary is nearly identical for 
galaxies divided into faint and bright samples (with a thresh- 
old of Mb = -18). The scatter between halo mass and galaxy 
luminosity smooths out the slight change in void radius with 
Mr seen in Figure |5] Thus for a luminosity threshold cut, the 
edge of the local void should be nearly the same for the two 
samples. 

The dashed curve i n Figure [5] indicate s the critical mass 
threshold below which iHoeft et al.l (l2006l) find that photoion- 
ization significantly reduces the ability of void halos to form 
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Fig. 4. — The void probability as a function of luminosity at three different 

separations relative to the mean intergalactic separation, Ig = ng . Points 
with errors are results from SDSS DR6. Curves represent results from the 
HOD model from the simulations Hsted in Table 2: L96 (dotted curves), L192 
(solid curves), L384 (dashed curves). Data points for the Mr = -22 sample are 
not shown for r > 1.3/« because the VPF is not measured at those distances. 



Fig. 5. — The maximum halo mass within voids as a function of dis- 
tance from the void center The solid curve represents mean of the largest 
three voids in the 96 /r' Mpc simulation. The points represent M,ni„ as a 
function of magnitude, placed along the curve to demonstrate how large the 
void would be at each luminosity. The dashed curve is the critical mass from 
Hoeft et al. (2006), below which star formation in void halos is significantly 
attenuated due to delayed formation histories of void halos. 



Stars. No halos atR<5 h~^ Mpc are above this line, suggest- 
ing that the inner half of deep voids will contain no galaxies at 
all, becoming truly dark. iGottlober et alj (l2003h calculated the 
mean void mass function for all halos within 80% of the void 
ra dius, finding good agre ement with the analytic predictions 
of ISheth & Tormeiil ([2002). However, they also find that the 
shape of the void mass function depends on position within 
the void itself. The mean void mass function in the L96 simu- 
lation (within 0.8/?void) is in good agreement with their results. 
It is important to consider scatter between mass and lumi- 
nosity for central galaxies. The points in Figure |5] represent 
Mmin, the mass at which the halo has only a 50% probability 
of containing a galaxy of that luminosity or brighter The error 
bars on each point indicate the mass at which that probability 
is 10% and if we double the scatter between mass and lumi- 
nosity to 0.3. Random fluctuations, or larger physical scat- 
ter between mass and luminosity, may place a galaxy slightly 
further toward the center of the void (or further away), but not 
enough to smooth out the sharp transition between filament 
and void seen in the simulations. 

4. SUMMARY AND DISCUSSION 

The void phenomenon consists of two observational facts: 
that voids contain few, if any, low-luminosity galaxies, and 
that the few void objects tend to have similar properties to the 
overall galaxy population. The controversial aspect is whether 
these facts are at odds with the current cosmology. Although 
the depth of voids and homogeneity of void objects are strik- 
ing features of the cosmic web, they are readily explainable 
within the context of ACDM, combined with a straightfor- 
ward model to connect galaxies and dark matter at all lumi- 
nosities and mass scales. 

The simple proposition within our implementation of the 



halo occupation distribution is that galaxy properties are de- 
termined solely by the mass of the halo in which the galaxy 
resides, independent of the halo's larger scale environment^ 
Although this model must break down at some high level of 
precision or with the details of certain galaxy properties, con- 
necting luminosity to halo mass is a robust method that has 
passed all tests thus far (lAbbas & ShethI 120061: ISkibba et all 
120061; [Tinker et al.l l2007al as well as the tests presented in 
this paper). The mass-only approach to galaxy bias readily 
explains both aspects of the void phenomenon. First, void 
galaxies are a fair representation of the field population sim- 
ply by construction. A galaxy in a 10'" /;"' M© halo does not 
know if the halo sits in a void or a filament. Recent numerical 
results have demonstrated that the halo itself does retain infor- 
mation about its environment (the 'assembly bias'; see, e.g., 
iGao et al. 2005; Wechsler et al. 2006; Gao & White 2007). 
However, this effect does not propagate into the void galax - 
ies through either luminosity or color (Tinker et al.ll2007"ah . 
iPatiri et al] (l2006bh find that void galaxy properties, even star 
formation rate or morphology, follow the same distribution 
as the field. Th i s is echoed in the semi-analytic results of 
ICroton & Farraii (2008). Therefore, although halos 'remem- 
ber' their environment, the physics of galaxy formation intro- 
duces an intrinsic scatter between halo formation and galaxy 
properties that washes this signal out. 

Second, the lack of dwarf galaxies in voids is a consequence 
of mapping galaxy luminosity to halo mass in such a way as to 
preserve the observed abundance of galaxies at each luminos- 
ity. The halo mass function scales as M~^, while the luminos- 

^ We note that it is possible to create an HOD with depends on environ- 
m ent. The model would then spe cify P(A^|M, <5) rather t han P(N\M), as done 
in lTinker et aj] <200ai2007if) and lWechslerltaJ] UM^ . 
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ity function scales as roughly L~' . This difference in logarith- 
mic slopes implies that the mass to light ratio M /L of galax- 
ies increases as galaxies become fainter Below L,, objects 
with substantially different luminosities will occupy halos of 
roughly comparable mass and will trace out the same large- 
scale structure. In order to explain the void phenomenon, 
we require the physical mechanism that causes the increase 
in M /L for faint objects (whatever it may be) to apply to all 
low-mass halos, regardless of environment. 

Connecting the luminosity function to the halo mass func- 
tion naturally predicts that the Tully-F isher relation becomes 
flatter at low luminosities. iGeha et aH (12006) confirmed this 
prediction by measuring the rotation velocities of galaxies as 
faint as M^ « -13. Our model places these galaxies in ha- 
los with maximum circular velocities of ^ 45 km s~' , in good 
agreement with the observed rotation speeds of 35-60 km s~' . 

As stated above, parameterizing the HOD as a function 
only of halo mass is only an approximation of the underly- 
ing physics (although a surprisingly good approximation). In 
the context of the void phenomenon, it should be noted that 
any influence of environment is mostly likely to reduce galaxy 
formation efficiency in underdense regi ons. This is the ef - 
fect of the photoionization arguments in iHoeft et aP (l2006|). 
and is also seen in the semi-analytic model of ICroton et al.l 



( |2007|) . A reduction in formation efficiency in voids makes 
^min higher at fixed luminosity. There is no need to in- 
voke exotic new physics, such as modified gravity models, to 
lower the formation efficiency of halos in underdense regions 
to match the data. Several studies have demonstrated that a 
long-range scalar potential acting on the dark matter can ac- 
celerate the evacuation of matter from voids during the growth 
of structure (iFarrar & Peeblesll2004t iGubser & Peebiell2004t 
iNusser et al.l 120051) . Although comparison of these models 
to observational data is difficult analytically, numerical sim- 
ulations of these a lternative gravity models, like those by 
iNusser et al.l (l2005h . can be used in the same manner as the 
simulations here. The excellent agreement between our HOD 
model and the observations may place upper limits on such 
models, since models that produce larger voids at any lumi- 
nosity would come into conflict with measured void statistics. 
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